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Abstract 
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' After a brief introduction to standard cosmology and the dark 

C*~> ■ matter problem in the Universe, we consider a self-gravitating nonin- 

, teracting fermion gas at nonzero temperature as a model for the dark 

| matter halo of the Galaxy. This fermion gas model is then shown 

(N) ■ to imply the existence of a supermassive compact dark object at the 

Galactic center. 
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Ph! 1 Introduction 

43 

At some stage of the evolution of the Universe, primordial density fluctuations 
must have become gravitationally unstable forming dense clumps of dark 
matter (DM) that have survived until today in the form of galactic halos. In 
the recent past, galactic halos were successfully modeled as a self-gravitating 
isothermal gas of particles of arbitrary mass, the density of which scales 
asymptotically as r -2 , yielding flat rotation curves ]]J. The aim of this paper 
is to describe the halo of our Galaxy in terms of a self-gravitating fermion 
gas in hydrostatic and thermal equilibrium at finite temperature. 

Self-gravitating weakly interacting fermionic matter has been exploited 
in a wide range of astrophysical phenomena. Originally, self-gravitating de- 
generate neutrino stars were suggested as a model for quasars 0, and later 
neutrino matter was used as a model for dark matter in galactic halos and 
clusters of galaxies, with a neutrino mass in the ~ eV range ||. Recently, 
degenerate superstars, composed of weakly interacting fermions in the ~ 10 
keV range, were suggested [Q, |5], [], [7], |8| as an alternative to the supermassive 
black holes that are believed to exist at the centers of galaxies. It was shown 
[|] that such degenerate fermion stars could explain the whole range of super- 
massive compact dark objects which have been observed so far, with masses 
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ranging from 10 6 to 3 x 1O 9 M , merely assuming that a weakly interacting 
quasistable fermion of mass m ~ 15 keV exists in Nature. Most recently, it 
has been pointed out that a weakly interacting dark matter particle in the 
mass range 1 % m/keV % 5 could solve the problem of the excessive structure 
generated on subgalactic scales in iV-body and hydrodynamical simulations 
of structure formation in this Universe ||. 

Of course, it is well known that the interval 1-15 keV lies squarely in 
the cosmologically forbidden mass range for stable active neutrinos v [TO]. 
The existence of such a massive active neutrino is also disfavored by the 



Super-Kamiokande data [|TTJ. However, for an initial lepton asymmetry of 
~ 10~ 3 , a sterile neutrino v s of mass m s ~ 10 keV may be resonantly pro- 
duced in the early Universe with near closure density, i.e., Q ~ 1 Hl2f . The 
resulting energy spectrum is not thermal but rather cut off so that it approx- 
imates a degenerate Fermi gas. In this mass range, sterile neutrinos are also 
constrained by astrophysical bounds on the radiative v s — > uj decay [OJ. 
However, the allowed parameter space includes m s ~ 15 keV, contributing 
n d ~ 0.3 to the critical density, as favored by the BOOMERANG data JTJ. 
As an alternative possibility, the ~ 15 keV sterile neutrino could be replaced 
by the axino [15| or the gravitino |16, [17| in soft supersymmetry breaking 
scenarios. 

As the supermassive compact dark objects at the galactic centers are 
well described by a degenerate gas of fermions, it is tempting to explore 
the possibility that one could describe both the supermassive compact dark 
objects and their galactic halos in a unified way in terms of a fermion gas at 
finite temperature. We will show that this is indeed the case, and that the 
observed dark matter distribution in the Galactic halo is consistent with the 
existence of a supermassive compact dark object at the center of the Galaxy 
which has about the right mass and size, and is in thermal and hydrostatic 
equilibrium with the halo. 



2 Standard cosmology 



Standard cosmology provides a successful description of the evolution of the 
Universe from a fraction of a second after the creation until today. A short 
review on the standard model of cosmology is given in | 18[| . For our purpose, 
it is sufficient to state the basic underlying principles. Standard cosmology 
is based on the following three theoretical assumptions: 

1. Cosmological principle. The cosmological principle asserts that the 
Universe is homogeneous and isotropic on large scales. The most general met- 
ric satisfying the cosmological principle is the Friedmann-Robertson- Walker 
metric ||TU| 



ds 2 = dt 2 -a(tf 



dr 2 



kr 2 



r 2 (d6 2 + sin dd<p 2 ) 



(1) 



where the curvature constant k takes on the values 1, 0, or -1, for a closed, 
flat, or open universe, respectively. The time-dependent quantity a{t) is the 
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scale factor of the expansion conveniently normalized to unity at present 
time, i.e., a(to) = 1. In other words, a is the radius of the Universe measured 
in units of its current radius. 

2. General relativity. Gravity is described by Einstein's general the- 
ory of relativity governed by the equivalence principle and Einstein's field 
equations. 

3. Perfect fluid. Matter is approximated by a homogeneous perfect 
fluid. The energy-momentum tensor then takes a simple form 

T^ v = (p + p)u,,u u - , (2) 

where p and p are the density and the pressure of the fluid, respectively. 

With these assumptions, the set of Einstein's equations reduces to the 
Friedmann-Robertson- Walker (FRW) equations 

- = ±-4n(p + 3p), (4) 
a 3 

where the natural system of units h = c = G = lis assumed. The first FRW 
equation describes the expansion of the Universe. The quantity H(t) is the 
Hubble "constant" and A is the cosmo logical constant. We define the critical 
density as p cr = 3ifo/87r an d the ratio of the density to the critical density 
is denoted by Q = p/p cr . The precise present value of the Hubble constant is 
not known, but the widely accepted value is H = 100h kms _1 Mpc _1 , with 
a dimensionless parameter h between 0.4 and 1. Dividing by H(t ) 2 , Eq. 
(|3[) at t = to can be conveniently written as a sum rule, i.e., 

n - k/(H a Q ) 2 + Q A = 1, (5) 

where Q A = A/(3Hq) is the vacuum energy contribution to the critical den- 
sity today. Observational evidence favors a flat universe today, i.e., k = 0. 
Thus, Eq. (H) becomes 

n + n A = i. (6) 

The second FRW equation describes the acceleration of the expansion. 
The expansion will accelerate or decelerate, depending on whether the vac- 
uum energy dominates the matter or vice versa. However, even if the cosmo- 
logical term vanishes, the expansion could accelerate if the dominant com- 
ponent of the DM obeys a peculiar equation of state such that the pressure 
is negative and higher than one third of the density. One popular example 
is the scalar field model called quintessence |[20|| . Another scenario is based 
on a fluid obeying the Chaplygin gas equation of state p oc — 1/p |f2T], E3|, 
which has been intensively investigated for its solubility in 1+1- dimensional 
space-time, for its supersymmetric extension and connection to <i-branes [E3 



The observational evidence in support of standard cosmology may be 



summarized in four empirical pillars ||24|| : Hubble's law, cosmic background 



radiation (CBR), anisotropy of CBR, and abundance patterns of light ele- 
ments. However, despite an overwhelming observational support, a number 
of problems remain unsolved: 
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• What caused the Big Bang and the expansion? Did the Universe begin 
with more than 3+1 dimensions? 

• Why is the cosmological constant A about 50-120 orders of magnitude 
smaller than the value expected from quantum field theory? 

• What caused the initial baryon-antibaryon asymmetry that led to the 
absence of antimatter today? 

• Why is the Universe so smooth on large scales as evidenced by CBR? 

• What caused the primordial density fluctuations that provided the 
seeds for structure formation? 

• What does nonbaryonic DM consist of? 

A solution to these problems most probably goes beyond the standard model 
of cosmology and certainly beyond the standard model of particle physics. 

3 Dark matter 

Here, we briefly discuss the DM problem and possible DM candidates. A 
more detailed analysis may be found in a number of recent review articles 

DM has to be introduced because of the following facts: 

• Astronomical observations, such as the flattness of the rotation curves 
of spiral galaxies and the peculiar motion of galaxies within clusters, 
strongly indicate that 

no _Anatt> 03 (?) 
Per 

• Consistency with the Big Bang nucleosynthesis implies for baryonic 
matter 

0.008/i 2 <n B < 0.024/i^ 2 . (8) 

• Astronomical observations yield a small relative density of luminous 
matter 

n lum = 0.003/io (9) 
From these facts we conclude that 

• About 95% of matter is dark. 

• About 60%-90% of DM is nonbaryonic. 

• At least 75% of baryonic matter is dark. 



4 



Whereas baryonic DM is most likely in the form of relatively standard as- 
trophysical objects, e.g., cold hydrogen clouds or compact objects such as 
neutron stars, brown dwarfs, MACHOs, and even black holes, the nature of 
nonbaryonic DM is unknown and still a subject of speculations. Candidates 
within the standard model are practically excluded and those beyond the 
standard model have not yet been detected in particle physics experiments. 
Nevertheless, cosmological and astrophysical observations tell us what these, 
yet undetected, particles could or could not be. 

The different DM scenarios are conveniently classified as hot, warm, and 
cold DM [ 26fl , depending on the thermal velocities of DM particles in the 



early Universe. 

Hot DM refers to low-mass neutral particles that are still relativistic when 
galaxy-size masses (~ 1O 12 M ) are first encompassed within the horizon. 
Hence, fluctuations on galaxy scales are wiped out by the "free streaming" of 
the dark matter. Standard examples of hot DM are neutrinos and majorons. 
They are still in thermal equilibrium after the QCD deconfinement transition, 
which took place at Tqcd — 150 MeV. Hot DM particles have a cosmological 
number density comparable with that of microwave background photons, 
which implies an upper bound to their mass of a few tens of eV. 

Warm DM particles are just becoming nonrelativistic when galaxy-size 
masses enter the horizon. Warm DM particles interact much more weakly 
than neutrinos. They decouple (i.e., their mean free path first exceeds the 
horizon size) at T ^> Tq CD . As a consequence, their number is expected to be 
roughly an order of magnitude lower and their mass an order of magnitude 
larger, than hot DM particles. Examples of warm DM are ~ keV sterile 
neutrinos, axinos or gravitinos in soft supersymmetry breaking scenarios 



16 , 17 . There has been renewed interest in the standard model neutrino as 



a candidate for warm DM 

Cold DM particles are already nonrelativistic when even globular cluster 
masses (~ 1O 6 M ) enter the horizon. Hence, their free streaming is of no 
cosmological importance. In other words, all cosmologically relevant fluctu- 
ations survive in a universe dominated by cold DM. The two main particle 
candidates for cold dark matter are the lowest supersymmetric weakly inter- 
acting massive particles (WIMPs) and the axion. 

One of the central issues in dark-matter modeling is the problem of struc- 
ture formation on subgalactic scales. The combination of cold DM and a 
small cosmological constant (ACDM) seems to be in good agreement with 
many observational constraints. However, N-body and hydrodynamical sim- 
ulations of galaxy formation evidence that ACDM overpredicts structure on 
small scales 0. In addition to that, high-resolution simulations generally 
find a dark matter profile with a central cusp p oc r~ 1,5 for galactic halos 
which seems to contradict the observations. 



Clustering on small scales could be suppressed by an upper limit to the 
phase-space density of DM particles owing either to degeneracy pressure if 
they are fermions or to a repulsive interaction if they are bosons. The fermion 
mass would have to lie in the range 0.1 m 10 keV. A specific scenario 
invoking keV mass fermions is the cool- dark-matter proposal [13j for which 
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candidates exist in shadow-world models and the axino. 



4 Galactic halo 

We now discuss the properties of the halo of our Galaxy assuming that it 
consists of a self-gravitating gas of keV mass fermions in hydrostatic and ther- 
mal equilibrium at finite temperature. The Milky Way Galaxy consists of 
five major components which are nested within each other f30fl . A spheroidal 
halo with modest concentations of stars and about 170 globular clusters ex- 
tends out to a radius of perhaps 200 kpc. Within a radius of ~ 25 kpc, the 
halo contains stars and open clusters that are concentrated into two essen- 
tially coplanar disks: the thin disk and the thick disk. In their innermost 
part the disks merge with a spheroidal bulge, the central concentration of 
luminous matter in the Galaxy. Finally, a compact dark object with a mass 
of M c ~ 2.6 x 10 6 M Q is located in the vicinity of the enigmatic radio source 
Sgr A* at the Galactic center ||31|| , within a radius of 18 mpc. 

Here, we demonstrate that by extending the Thomas-Fermi theory to 
nonzero temperature it is possible to explain, within the same model, both 
the Galactic halo and the compact dark object at the Galactic center. Ex- 
tending the Thomas- Fermi theory to finite temperature |32], [33], |34j], it has 



been shown that, at some critical temperature T c , weakly interacting mas- 
sive fermionic matter with a total mass below the Oppenheimer-Volkoff limit 



[ |35|1 undergoes a first-order gravitational phase transition from a diffuse to a 
clustered state, i.e., a nearly degenerate fermion star. However, during this 
first-order phase transition a large amount of latent heat must be released 
in order to substantially decrease the entropy of the initial diffuse config- 
uration. In the absence of a mechanism which would make such a release 
possible, the system would remain in a thermodynamic quasistable super- 
cooled state close to the point of gravothermal collapse. The Fermi gas will 
be caught in the supercooled state even if the total mass of the gas exceeds 
the Oppenheimer-Volkoff limit as a stable condensed state does not exist in 
this case. 

The formation of a quasistable supercooled state may be understood as 
a process similar to that of violent relaxation, which was introduced to de- 



scribe rapid virialization of stars of different mass in globular clusters [36, ^7 
Through the gravitational collapse of an overdense fluctuation at about one 
Gyr after the Big Bang, part of gravitational energy transforms into the 
kinetic energy of random motion of small-scale density fluctuations. The re- 
sulting virialized cloud will thus be well approximated by a gravitationally 
stable thermalized halo. In order to estimate the mass-to-temperature ratio, 
we assume that an overdense cloud of mass M stops expanding at the time t m , 
reaching its maximal radius R m and the average density p m = 3M/(47ri?^). 
The total energy per particle is just the gravitational energy 

^ 3 1 
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P(r) = —■ (12) 



From the spherical model of nonlinear collapse it follows 

9tt 2 _ 9tt 2 

Pm = — p(t m ) = —VL d p cr (l + Z m ) 3 , (11) 

lb lb 

where p(t m ) is the background density at the time t m or the cosmological 
redshift z m . We approximate the virialized cloud by a singular isothermal 
sphere [^] of the mass of the Galaxy M and radius R. The singular isother- 
mal sphere is characterized by a constant circular velocity O = (2T/m) 1//2 
and the density profile 

_0 2 

47rr 2 

Its total energy per particle is the sum of gravitational and thermal energies, 

£ = -I^ = -V (13) 
4 R 4 v ; 

Combining Eqs. Qig), QTJ, and (0), we find 

6 2 = ^(b^p cr M 2 ) 1 / 3 (l + 2 m ). (14) 

Taking fi d = 0.3, M = 2 x 1O 12 M , z m = 4, and h = 0.65, we find ~ 
220 km s -1 , which corresponds to the mass-temperature ratio m/T ~ 4 x 10 6 . 
Next, we briefly discuss the general-relativistic Thomas- Fermi theory |34] 



for a self-gravitating gas, consisting of iV fermions of mass m in equilibrium 
at a temperature T, enclosed in a sphere of radius i?. We denote by p, 
p, and n the pressure, energy density, and particle number density of the 
gas, respectively. The metric generated by the mass distribution is static, 
spherically symmetric, and asymptotically flat, i.e., 

ds 2 = i 2 dt 2 - (1 - 2M/r)- 1 dr 2 - r 2 [dQ 2 + sin#c/0 2 ). (15) 

For numerical convenience, we introduce the parameter 

a = | (16) 

and the substitution 

£=^ + l)- 1/2 , (17) 
m 

where \x is the chemical potential associated with the conserved particle num- 
ber N. Using this, the equation of state for a self-gravitating ideal gas may 
be represented in a parametric form ||39|| : 



1 Z" 00 y 2 
U= ^Jo dV l + exp{[(|/ 2 + iy/ 2 /{cp+ l) 1 / 2 -!]«}' MlN 



7T 2 



y2/ 2 , 1 U/2 

dv - W ' (19) 

y l + exp{[(y 2 + iy/ 2 /{p + 1)V2 -l]a}' 1 J 

= J_ r . i/ 4 (i/ 2 + iy 1/2 (9(]] 

P 37r 2 io V l + exp{[{y 2 + iy/y{<p + 1)V2 -l]a}' 1 J 
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We have chosen appropriate length and mass scales a and 6, respectively, 
such that 

a = h = \h^ (21) 



9 ' 

g m z 



or, restoring H, c, and G, we have 



f2hM Fl 1 nrmn , M0 /2 /!5keV\ 



2 



1.0798 x 10 1U J- km, (22) 




g cm? 



^Mk = 0.7251 x 10 10 J- I ) M . (23) 

;/ ///- ' 



Here, Mpj = J%c/G denotes the Planck mass and g the combined spin- 
degeneracy factor of neutral fermions and antifermions, i.e., g=2 or 4 for 
Majorana or Dirac fermions, respectively. In this way, the fermion mass, 
the degeneracy factor, and the chemical potential are eliminated from the 
equation of state. 

Einstein's field equations for the metric fllS]) are given by 

i = - 2 ^ +1) 7^M)' (24) 



dr 

To these two equations we add 

dN 



dM n , , 

Anr 2 p. (25) 



47rr 2 (l - 2M/r)- 1/2 n, (26) 
dr 

imposing the particle-number constraint as a condition at the boundary 

AT{R) = N. (27) 

Equations (H)-(H) should be integrated using the boundary conditions at 
the origin: 

^(0) =<p > -1; A^(0) = 0; 7V(0) = 0. (28) 

It is useful to introduce the degeneracy parameter rj = ap/2, which, in the 
Newtonian limit, approaches r/ nr = (yU nr — V)/T. Here, we have introduced 
the nonrelativistic chemical potential /i nr = fi — m, with fi nT -C m, and the 
approximation £ = e y l m ~ 1 + V/m, with V being the Newtonian potential. 
As (p is monotonously decreasing with increasing r, the strongest degeneracy 
is obtained at the center with 770 = aipo/2. The parameter 770, uniquely 
related to the central density and pressure, will eventually be fixed by the 
requirement (p7| ). For r > R, the function ip yields the usual empty-space 
Schwarzschild solution 

a 2 ( 2M\- 1 , . 

^=£3(1 -1, 29 

m 2 \ r J 
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Figure 1: Energy (shifted by 12 x 10 8 m) versus temperature (shifted by 
-24 x 10~ 8 m), both in units of lO" 10 ™, for fixed JV = 2x 10 12 M o /m 

with 

r R 

M = M{R)= drA-nr 2 p(r). (30) 
Jo 

Given the temperature T, the set of self-consistency equations (|i"8|)-(^6|), with 
the boundary conditions ( p7|) -(|3"0D defines the general-relativistic Thomas- 
Fermi equation. 

The numerical procedure is now straightforward. For a fixed, arbitrarily 
chosen a, we first integrate Eqs. ([24]) and (p5|) numerically on the interval 
[0, R] and find solutions for various central values t]q. Integrating fl26f) simul- 
taneously, we obtain J\f(R) as a function of r] . We then select the value of 
?7o for which J\f(R) = N. The chemical potential \x corresponding to this 
particular solution is given by (p9|). If we now eliminate \i using (|16[), we 
finally get the parametric dependence on temperature through a. 

The quantities N, T, and, R are free parameters in our model and their 
range and choice are dictated by physics. At T = the number of fermions 
N is restricted by the OV limit iVov = 2.89 x 10 9 (15 keV /m) 2 M Q /m. 
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However, at nonzero temperature, stable solutions exist with iV > Nov, 
depending on temperature and radius. In the following, N is required to be 
of the order 2 x !O 12 M /m, so that for any m, the total mass is close to 



the estimated mass of the halo |]4"0 |. As we have demonstrated, the expected 
particle mass-temperature ratio of the halo is given by a ~ m/T = 4 x 10 4 . 
The halo radius R is in principle unlimited; in practice, however, it should not 
exceed half the average intergalactic distance. It is known that an isothermal 
configuration has no natural boundary, in contrast to the degenerate case of 
zero temperature, where for given iV (up to the OV limit) the radius R is 
naturally fixed by the condition of vanishing pressure and density. At nonzero 
temperature, with R being unbounded, our gas would occupy the entire 
space, and fixing iV would make p and p vanish everywhere. Conversely, if 
we do not fix N and integrate the equations on the interval [0, oo), both M 
and iV will diverge at infinity for T > 0. Thus, one is forced to introduce a 



cutoff. In an isothermal model of a similar kind the cutoff was set at 
the radius R, where the energy density was by about six orders of magnitude 
smaller than the central value. Our choice of R = 200 kpc is based on the 
estimated size of the Galactic halo. 

The only remaining free parameters of our model are the fermion mass m 
and the degeneracy factor g, which always appear in the combination m 4 g. 
We fix these parameters at m = 15 keV and g = 2, and justify this choice a 
posteriori. 

We now present the results of the calculations for fixed particle number 
and temperatures near the point of gravothermal collapse. In Fig. [I] the 
energy per particle defined as E = M/N — m is plotted as a function of 
temperature for fixed N = 2 x lO 12 M /m. The plot looks very much like 
that of a canonical Maxwell-Boltzmann ensemble J37j, with one important 
difference: in the Maxwell-Boltzmann case, the curve continues to spiral in- 
wards ad infinitum approaching the point of the singular isothermal sphere, 
that is characterized by an infinite central density. In the Fermi-Dirac case, 
the spiral consists of two, almost identical curves. The inwards winding of 
the spiral begins for some negative central degeneracy and stops at the point 
T = 2.3923 x 10~ 7 m, E = —1.1964 x 10 _7 m, where rjo becomes zero. This 
part of the curve, which basically depicts the behavior of a nondegenerate 
gas, we call the Maxwell-Boltzmann branch. By increasing the central de- 
generacy parameter further to positive values, the spiral begins to unwind 
outwards very close to the inwards winding curve. The outwards winding 
curve will eventually depart from the Maxwell-Boltzmann branch for tem- 
peratures T ^ 10 _3 m. Further increase of the central degeneracy parameter 
brings us to a region where general-relativistic effects become important. The 
curve will exhibit another spiral for temperatures and energies of the order 
of a few 10~ 3 m approaching the limiting temperature = 2.4 x 10~ 3 m 
and energy = 3.6 x 10~ 3 m, with both the central degeneracy parameter 
and the central density approaching infinite values. It is remarkable that 
gravitationally stable configurations with arbitrary large central degeneracy 
parameters exist at finite temperature even though the total mass exceeds 
the OV limit by several orders of magnitude. 
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Figure 2: Number of particles versus central degeneracy parameter for 
m/T = 4 x 10 6 (solid), 3.5 x 10 6 (short dashs), 4.5 x 10 6 (long dashs), and 
5 x 10 6 (dot-dashed line). 



The results of the numerical integration of Eqs. (p4|) and (pS"l), without 
restricting N, are presented in Fig. |2|, where we plot the particle number N 
as a function of the central degeneracy parameter 770 for several values of 
a close to 4 x 10 6 . For fixed N, there is a range of a, where the Thomas- 
Fermi equation has multiple solutions. For example, for N = 2 x 10 12 and 
a = 4 x 10 6 six solutions are found, which are denoted by (1), (2), (3), 
(3'), (2'), and (1'), corresponding to the values rj Q = -30.528, -25.354, - 
22.390, 29.284, 33.380, and 40.479, respectively. In Figs. | and | we plot the 
corresponding density profiles and enclosed masses, respectively. For negative 
central value 770, for which the degeneracy parameter is negative everywhere, 
the system behaves basically as a Maxwell-Boltzmann isothermal sphere. 
Positive values of the central degeneracy parameter 770 are characterized by 
a pronounced central core of mass of about 2.5 x 1O 6 M within a radius of 
about 20 mpc. The presence of this core is obviously due to the degeneracy 
pressure of the Fermi-Dirac statistics. A similar structure was obtained in 
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Figure 3: The density profile of the halo for a central degeneracy parameter 
t] = (dotted line) and for the six ^-values discussed in the text. Configu- 
rations with negative r] ((l)-(3)) are depicted by the dashed and those with 
positive r] Q ((l')-(3')) by the solid line. 



collisionless stellar systems modeled as a nonrelativistic Fermi gas [42 



Figs. |3| and |^ show two important features. First, a galactic halo at a given 
temperature T may or may not have a central core, depending on whether 
the central degeneracy parameter t] is positive or negative. As the potential 
is nearly harmonic up to about 1 to 10 kpc for negative r] , this may favor the 
formation of a barred galaxy. Second, the closer to zero rjo is, the smaller the 
radius is at which the r~ 2 asymptotic behavior of the density begins. The 
flattening of the Galactic rotation curve begins in the range 1 ^ r/kpc 1z 10, 
hence the solution (3') most likely describes the Galaxy's halo. This may be 
verified by calculating the rotational curves in our model. We know already 
from our estimate (|i~4]) that our model yields the correct asymptotic circular 
velocity of 220 km/s. In order to make a more realistic comparison with 
the observed Galactic rotation curve, we must include two additional matter 
components: the bulge and the disk. The bulge is modeled as a spherically 
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Figure 4: Mass of the halo Mh(r) enclosed within a radius r for various 
central degeneracy parameters rj as in Fig. [3|. 



symmetric matter distribution of the form [[43 1 



e -/is „oq e ~hsu 

Ph{s) = ^l du WTW^W J (31) 

where s = (r/ro) 1 ^ 4 , ro is the effective radius of the bulge and h is a parame- 
ter. We adopt ro = 2.67 kpc and h yielding the bulge mass Mb = 1.5x 10 10 M Q 
[ [yj] . In Fig. H the mass of halo and bulge enclosed within a given radius is 
plotted for various 770. Here, the gravitational backreaction of the bulge on 
the fermionic halo has been taken into account. The data points, indicated 
by squares, are the mass M c = 2.6 x 1O 6 M within 18 mpc, estimated from 

h0.2 v 1n ll 



the motion of the stars near Sgr A* ||45|| , and the mass M 50 = 5.4_ 3 ; 6 x 10 



within 50 kpc, estimated from the motions of satellite galaxies and globular 
clusters [fiO|l . Variation of the central degeneracy parameter r] between 24 



and 32 does not change the essential halo features. 

In Fig. |6] we plot the circular velocity components of the halo, the bulge, 
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Figure 5: Enclosed mass of halo plus bulge versus radius for 7] 
(dashed), 28 (solid), and 32 (dot-dashed line). 



24 



and the disk. The contribution of the disk is modeled as 



d (r o ) 



1.97(r/r c 



,1.22 



(r/r ) 2 + 0.78 2 ] 1 - 43 



(32) 



where we have taken r Q = 13.5 kpc and 0d = 100 km/s. Here it is assumed 
for simplicity that the disk does not influence the mass distribution of the 
bulge and the halo. Choosing the central degeneracy r/o = 28 for the halo, 



the data by Merrifield and Oiling |57j are reasonably well fitted. 

We now turn to the discussion of our choice of the fermion mass m = 15 
keV for the degeneracy factor g = 2. To that end, we investigate how the 
mass of the central object, i.e., the mass M c within 18 mpc, depends on m 
in the interval 5 to 25 keV, for various rjo. We find that m ~ 15 keV always 
gives the maximal value of M c ranging between 1.7 and 2.3 xlO 6 M for 770 
between 20 and 28. Hence, with m ~ 15 keV we get the value closest to the 
mass of the central object M c estimated from the motion of the stars near 
SgrA* 



45 
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Figure 6: Fit to the rotation curve of the Galaxy. The data points are from 



47j for R = 8.5 kpc and 6 = 220 km/s. 



5 Conclusions 

In summary, using the Thomas-Fermi theory, we have shown that a weakly 
interacting fermionic gas at finite temperature yields a mass distribution 
that successfully describes both the center and the halo of the Galaxy. For a 
fermion mass m ~ 15 keV, a reasonable fit to the rotation curve is achieved 
with the temperature T = 3.75 meV and the central degeneracy parameter 
?7o = 28. With the same parameters, we obtain the mass M50 = 5.04 x 10 U M & 
and M 2 oo — 2.04 x 1O 12 M within 50 and 200 kpc, respectively. These values 
agree quite well with the mass estimates based on the motions of satellite 
galaxies and globular clusters [|D|. Moreover, the mass of M c ~ 2.27 x 



1O 6 M , enclosed within 18 mpc, agrees reasonably well with the observations 
of motion of stars near the compact dark object at the center of the Galaxy. 
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